How to Falsify the GR+ACDM Model 
with Galaxy Redshift Surveys 
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A wide range of models describing modifications to General Relativity have been proposed, but no 
fundamental parameter set exists to describe them. Similarly, no fundamental theory exists for dark 
energy to parameterize its potential deviation from a cosmological constant. This motivates a model- 
independent search for deviations from the concordance GR+ACDM cosmological model in large 
galaxy redshift surveys. We describe two model-independent tests of the growth of cosmological 
structure, in the form of quantities that must equal one if GR+ACDM is correct. The first, e, 
was introduced previously as a scale-independent consistency check between the expansion history 
and structure growth. The second, v, is introduced here as a test of scale-dependence in the 
linear evolution of matter density perturbations. We show that the ongoing and near-future galaxy 
redshift surveys WiggleZ, BOSS, and HETDEX will constrain these quantities at the 5 — 10% level, 
representing a stringent test of concordance cosmology at different redshifts. When redshift space 
distortions are used to probe the growth of cosmological structure, galaxies at higher redshift with 
lower bias are found to be most powerful in detecting the presence of deviations from the GR+ACDM 
model. However, because many dark energy or modified gravity models predict consistency with 
GR+ACDM at high redshift, it is desirable to apply this approach to surveys covering a wide range 
of redshifts and spatial scales. 

PACS numbers: 



I. INTRODUCTION 

Evidence for accelerated cosmic expansion poses 
a serious challenge to the standard paradigm of Big Bang 
cosmology, where gravity is described by General Rela- 
tivity (GR), and the Universe is filled by matter and ra- 
diation fields. Most proposed explanations invoke failure 
of GR on cosmological scales or introduce a new com- 
ponent with strongly negative equation of state in the 
cosmic energy budget. Without compelling theoretical 
reasons to support any particular theory, Occam's razor 
favors the model able to explain all of the cosmological 
observations with the minimum number of parameters. 
Currently, that model is "GR+ACDM" , characterized by 
two assumptions: gravity is correctly described by GR, 
and cosmic acceleration is caused by a cosmological con- 
stant (A). We refer to it as the concordance cosmological 
model. 

The rich literature in this field offers many explana- 
tions for the accelerated expansion alternative to the con- 
cordance model (for a list of references, see [3]). Accord- 
ing to which of the two key assumptions is broken, they 
fall in the category of Modifications of Gravity (MoG) 
or Dark Energy (DE). Instead, we consider the concor- 
dance model in its entirety as the reference framework 
and ask ourselves how we can most thoroughly test that 
all cosmological observations are in agreement with it. 

The motivation for this approach is manifold. First 
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of all, it is feasible. While interpretation of the data 
in the context of a given theory is necessarily model- 
dependent, the detection of a significant deviation from 
an expected behavior is incontrovertible. Moreover, the 
unlimited freedom intrinsic in more complicated MoG 
and DE models is reflected by degeneracies in cosmologi- 
cal observables, so that it is impossible to uniquely iden- 
tify the phenomenology of one or the other class, without 
making further theoretical assumptions {e.g., [4]). Sec- 
ond, our approach is relevant. If the concordance model 
were proven to be faulty, this would profoundly change 
our understanding of the Universe, whatever the exact 
explanation might be. Finally, this approach is sensible; 
practicality suggests that if the concordance model is fal- 
sified it would motivate a far more detailed exploration 
of the possible underlying reasons, a search then guided 
by knowledge of the nature, and location in space and 
time, of the observed deviation. 

The paper is organized as follows. The first step, dis- 
cussed in Sec. HH is to identify relations between cos- 
mological observables, which hold for the concordance 
model. This is similar in spirit to what was proposed 
in [5|, but here we further break down the concordance 
model behavior in order to obtain two independent re- 
lationships. This allows us to identify two parameters, 
named e and v, whose value can be exactly predicted 
for the GR+ACDM case. In Sec. OH] we write these 
parameters in terms of observable quantities, and we se- 
lect the methods for measuring the quantities entering 
the definition of e and v. We focus on the use of the 
power spectrum and bispectrum of galaxy redshift sur- 
veys, in conjunction with CMB measurements. The es- 
tablished framework is used in Sec. II VI to define a Figure 
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of Merit in the (e, v) plane for the task at hand. The 
computation of the Figure of Merit requires knowledge 
of the observational uncertainties in the measurement 
of e, v, and the other relevant cosmological parameters; 
these are determined using the Fisher Matrix technique. 
We compute the actual Figure of Merit for the surveys 
WiggleZ, BOSS and HETDEX. This approach allows a 
direct model-independent comparison of the capabilities 
of different experiments in answering the question "Is 
GR+ACDM the correct cosmological model?" We sum- 
marize our results and list concluding remarks in Sec. 

El 



II. THEORETICAL FRAMEWORK 

We assume scalar linear perturbations around a flat 
FRW background in the Newtonian gauge, 



ds 2 



(1 + 2^)dt 2 + a 2 (l - 2§)dx 2 



(1) 



and work in the quasi-static, linear approximation, which 
is valid for sub-horizon modes still in the linear regime 

e.g., !• 

The evolution of perturbations is described by the con- 
tinuity, Euler and Poisson equations {e.g., [7|): 



6' = -—v 

aH ' 



V m + V m = — — * 

aH 



/c 2 $ = -AirGa 2 p m £„ 



(2) 

(3) 
(4) 



where S m and v are the matter density perturbations 
and the divergence of the peculiar velocity, respectively. 
In the concordance model, the anisotropic stress van- 
ishes, and these equations can be combined to derive 
the equation of motion for the matter density pertur- 
bations, or equivalently for the growth factor D defined 
as D(k,a) — S m (k,a)/S m (k,a = 0): 



D" + (2 + %)D' 
H 



AttG 



p m D = 0; 



(5) 



where a is the scale factor, H = a/a is the Hubble func- 
tion, a prime denotes derivative with respect to In a, and 
p m is the average matter density. 

Two well-known features of the GR+ACDM model are 
immediately evident from the form of this equation (e.g. 
[a]). The solution for the linear growth factor can be 
determined exactly once the Hubble function is known; 
and, since none of the coefficients is a function of scale, 
the growth factor is scale- independent. These are two 
independent properties of the concordance model, and we 
propose to test them separately using the two parameters 



e(a) 
v(a) 



fV (a) /(«); 

f(k s ,a) 
f{kh,a)' 



(6) 
(7) 



In the above equations, Q rn is the matter contribution 
to the total density, / = , fc s ( ma ii) and fc L (ar g o) rep- 
resent two different scales where the logarithmic deriva- 
tive of the growth is measured, and 7(a) is the sim- 
ple function entering the well known fitting function 
(dlnD/dlno)fit - [rj m (a)] 7(a) . Throughout this paper, 
we use for the fit the value 7 = 0.55 [8j, which is accu- 
rate at the 0.05% level for GR+ACDM at the redshifts 
of interest. 

The definition of the parameter e is similar to the one 
used in Q but with one important difference: we as- 
sume here that the measurement of / in e is averaged 
over all the scales probed by each galaxy redshift sur- 
vey. Therefore, the parameter e is only sensitive to the 
consistency between the Hubble expansion (equivalent to 
^m(a)) and the growth as a function of time D(a), re- 
quired by the concordance model, while v is sensitive to 
a scale dependence in the logarithmic derivative of the 
growth. If many measurements of the growth as a func- 
tion of scale are available, v could be defined in multiple 
ways. We choose to maximize the signal to noise for the 
measurement of this parameter by using two bins in k 
for each survey, and choosing the ranges so that the two 
measurements have similar fractional uncertainties. 

Both parameters are exactly = 1 at all times in the 
GR+ACDM model, and this formalism allows for an easy 
detection of deviations from the concordance model be- 
havior, although it would not tell us immediately whether 
non-A Dark Energy or MoG is the culprit. 

A scale dependence in the logarithmic derivative of the 
growth factor (corresponding to D ^ 1) is challenging 
within the context of perfect fluid Dark Energy models 
0, but is predicted by models with varying sound speed, 
or anisotropic stress e.g., [IollTTl|. Limits on such quanti- 
ties are still weak [l2[ • Massive neutrinos may also induce 
a scale-dependent growth factor. Neutrinos are assumed 
massless in the concordance model. In the case of neutri- 
nos of non- negligible but known mass, the effect of such 
scale dependence can be removed and the formalism ap- 
plied. If the same galaxy survey is used to measure the 
logarithmic growth and the mass of neutrinos, we expect 
the Figure of Merit to be smaller. 



III. FROM THEORY TO OBSERVABLES 

In order to compute a Figure of Merit for the measure- 
ment of e, and v, we start by rewriting them in terms of 
observable quantities. We have: 

e(a) = Q-^f(a) = f" H{ f" B(a)b(a)- (8) 

v ( a ) = ^( fcs ' a ) = l 3 ( k s, a )b( a ) = P(k s ,a) , g s 
f(k L ,a) fi(k L .a)b(a) /3(k L ,a)' 

Here, we have written / as /3 • 6, where /3 is the linear 
redshift distortion parameter, and b is the bias of the 
target galaxies. We assume that any scale dependence 
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of the bias has been removed; see Sec. IIVI for further 
discussion. 

CMB experiments are an excellent probe of the combi- 
nation il mi oiJQ, whose value is known to within 5% from 
WMAP7 data [l3j . We assume that this measurement 
will be performed by the ongoing CMB satellite mission 
Planck [1J] with a relative precision of 1.5%. We repre- 
sent this constraint as a Gaussian prior in the likelihood 
function. 

The linear bias of galaxies can be measured in several 
ways; here we assume it has been measured via the 
bispectrum of galaxies, as suggested in [l5j]. That 
paper reports the precision for the bias measurement 
achievable by several reference surveys, together with 
a discussion on how the precision varies as a function 
of survey parameters (e.g., number density of galaxies, 
volume, fcmax)- We use simple scaling relations, de- 
rived from their discussion, in order to determine how 
well galaxy surveys can measure the bias, and again 
introduce this constraint in the likelihood function as 
a Gaussian prior. We assume that no information on 
the bias is derived from the power spectrum, since a 
correct implementation of a simultaneous measurement 
of the bias from bispectrum and power spectrum would 
require knowledge of the covariance between them. We 
note that as a result, we might be underestimating the 
potential of a given survey in measuring the bias. 
The remaining quantities appearing in the definition of e 
and v are the Hubble parameter at a given time, H(a), 
and the linear redshift distortion parameters measured 
on two different scales, f3(k s ,a) and P(kL,a). These 
functions can be measured by galaxy redshift surveys. 
It was shown in [16| that fitting the full 2-D galaxy 
power spectrum allows one to maximize the signal to 
noise in the measurement of H (a) and (3. We modify the 
public Fisher Matrix code available at [171 to account 
for the scale dependence of the linear redshift distortion 
parameter and to adopt the new parameterization in 
terms of e and v. 



primordial perturbations, and a v is the peculiar velocity 
disp ersion we use to model the Finger of God effect 
[Hf, following [H, 03. We use the measurement of 
the galaxy power spectrum in redshift space for the 
first six parameters, and assume no covariance between 
them and the measurements of the bias (obtained from 
bispectrum) and n m ^Hg (obtained from the CMB), as 
described above. 

For a galaxy redshift survey, the Fisher Matrix has the 
form (e.g., [20|) 



F, 



4irk 2 dk r 1 , <91nP„(fc,u) dlnP„(k,u) „ , 

-w(k,jJL). 



(2tt) 3 



dp- ., 

opi 



d Pj 



(10) 

Here \x is the cosine of the angle between the vector k 
and the line of sight, w(k,fi) is the weight function 



w(k,n) = - 



n g P g (k,n) 



l + n g P g (k,fj,) 



V. 



(11) 



accounting for the survey's volume and number density 
of galaxies n g , and /c max is the maximum wavenumber for 
which the formalism is valid |57[. We model the linear 
part of P g (k, fi) analytically as in [Ta |: 

P g (k,fi) = b 2 [P ss (k) + 2fri 2 P Sv (k) + p 2 ^P vv (k)]; (12) 

S refers to density and v to velocity perturbations. 

Many of the relevant derivatives appearing in the 
Fisher Matrix formula (|10l) are given in the appendix 
of [16|, and we do not report them here. The newly in- 
troduced derivatives with respect to (3 S and Pl can be 
obtained using the same formula as the derivative with 
respect to ft, but integrating in two separate k bins. If 
fc me d is the wavenumber chosen to separate the measure- 
ment of "small" and "large" scales, 



d\nP g (k) _ afu L ^lnP(fc) 



dp s 
d\nPg(k) 

d(3 L 



— 8(k, fc m od) " 



d/3 



(13) 



(l-0(fc,fc mod ))^^; (14) 



IV. FIGURE OF MERIT 

We now perform a Fisher Matrix calculation and use 
the results to generate a Figure of Merit in the plane (e, 
v). We assume that each survey provides a measurement 
of the relevant parameters at a single redshift (typically 
the mean redshift of the survey) , and therefore we drop 
the time dependence of the parameters. This allows one 
to directly compare the performance of different surveys, 
even when they target different redshift ranges. For no- 
tation convenience, we denote P(k s ) and P(k^) as f3 s and 
Pl respectively. 

As a first step, we compute the Fisher 
Matrix for the set of parameters pi = 
{P s ,P L ,\nH,\nA,n s ,b,a v ,n m . H^}; A and n s rep- 
resent the amplitude and scalar spectral index of the 



where 9(k\, k2) — if k\ < k2, and 1 otherwise. 

The second step is to transform to the final set of 
parameters qi = {e, v, hiH, hi A, n s , b, a v , f2 m .o-ffo}' an< l 
compute the corresponding Fisher Matrix, F a \>: 



F,- V^^F 



(15) 



For this purpose, it is useful to rewrite the parameters e 
and v explicitly in terms of /3 S , Pl, H, b and fl m fiHQ: 



Pi 



(16) 
(17) 



where we have used the fact that the parameter P ap- 
pearing in e is an average of the measurements over small 
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Galaxy 




V 




^mcd 


fcmax 


b 


A6/6 


Survey 




(h- 3 Gpc 3 ) 




(fl/Mpc) 


(h/Mpc) 






WigglcZ 


0.24 ■ 10 b 


1.0 


0.6 


0.1 


0.2 


1 


0.06 


BOSS 


1.5 ■ 10 b 


4.4 


0.45 


0.1 


0.18 


2 


0.03 


HETDEX 


0.8 ■ 10 B 


3.0 


2.7 


0.15 


0.4 


2 


0.015 



TABLE I: Specifics of the considered surveys. We display 
the total number of galaxies N &a ,\, the volume V, the mean 
redshift z m , the wavenumber fc mc d where we split the surveys, 
the maximum linear wavenumber fcmax, the galaxy bias b, and 
the percentage error in the bias. 



and large scales, a is the scale factor at the mean redshift 
of each survey, and 7 is the corresponding value of the 
fitting function introduced earlier. For the concordance 
model, p s = f3 L = (3. 

The non-null derivatives appearing in Eq. (|15[) are: 



1 1 1 1 1 1 



BOSS 

WiggleZ 

HETDEX 




de 

dlnH 
de 

d(£l m< oH ) 
de 

db 
de 

dv 
dp L 
dv 



dfi L = 2 jiW^, 
de b 
( q'fl* i~ 7 



7 



A+/3z 



(18) 
(19) 

(20) 

(21) 
(22) 
(23) 



With these derivatives in hand, we can now compute 
the Figure of Merit for three representative ga laxy sur- 
veys, WiggleZ (ongoi ng 1211 ) . BOSS (ongoing, [22j ) and 
HETDEX (planned, [23. 12411. t argeting respectively blue 
emission-line galaxies (e.g., [251] 1. Luminous Red Galaxies 
(e.g., [HI) and Lyman Alpha Emitters (e.g., [27J). The 
assumed specifics of the surveys are reported in Table Q] 
the value of Ab is the la amplitude of the Gaussian prior 
on the linear bias. For WiggleZ, we assume the final size 
of the survey as a reference. 

The marginalized errors in the (e,v) plane for the three 
surveys are shown in Figure CD While the three surveys 
are exploring different redshift windows, it is fair to plot 
the constraints on the same plane, because e and v are 
constructed to be exactly equal to one for the concor- 
dance model at all redshifts. The inverse of the area of 
the ellipse thus provides a Figure of Merit that allows 
a fair comparison of the strength of different surveys in 
detecting deviations from the concordance model. How- 
ever, constraints centered at different redshifts are inde- 
pendent from each other. 



FIG. 1: Expected errors in the (e,v) plane for the WiggleZ, 
BOSS and HETDEX surveys. The concentric ellipses show 
the 1 and 2 a contours. The values of the Figure of Merit, 
defined as n/A, where A is the area of the 1-a ellipse, are 204, 
93 and 262, respectively. 



Figure Q] shows that HETDEX will perform slightly 
better than WiggleZ in detecting deviations from the 
GR+ACDM model, and both will improve significantly 
on the result from BOSS. More specifically, the la uncer- 
tainty in the measurement of e, v is 6%, 4% for WiggleZ, 
10%, 9% for BOSS, and 5%,4% for HETDEX. The cor- 
responding values of the Figure of Merit, defined as tt/A, 
are 204, 93 and 262 respectively. 

The results mainly depend on the following factors. 
First, the measurement of the linear bias obtained via 
the bispectrum is noticeably better at higher redshift, as 
evident from Figs. 5 and 6 in |l5j |. Higher- redshift sur- 
veys are also favored by the larger number of available 
modes in the linear regime. Similarly, the signature of 
d\n D /din a (equivalent to j3 for galaxies of same bias) 
is larger at z ~ 2 than z ~ 0.5, because the growth 
is more rapid during matter domination than at later 
times, when the dark energy slows down the growth of 
the perturbations. These three factors reduce the un- 
certainties in the measurement of e and v, in the case of 
HETDEX vs. BOSS. Finally, for a similar rate of growth, 
the linear redshift distortion signal is largest for galaxies 
of low bias (such as the ones targeted by WiggleZ). This 
not only positively influences the measurement of j3 s and 
Pl, but also alters the correlation between e and v. For 
the particular configuration of WiggleZ, e and v are less 
correlated than for the other two surveys. This effect 
is slightly mitigated by the worse sensitivity of WiggleZ 
to a v and /?, resulting from the smaller volume and bias 
factor, but the resulting Figure of Merit is still excep- 
tionally promising. One possible caveat about the use of 
the lowest-bias galaxies comes from the N-body Simula- 
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tions performed by the authors of [28| . We have assumed 
here that the scale dependence of the bias is either ab- 
sent, or it is possible to correct for it. Their work shows 
that the bias is reasonably independent of scale in the 
most massive halos, where high-bias galaxies reside, but 
not in low-mass halos, where low-redshift Emission Line 
Galaxies are. If accounting for such scale dependence 
turns out to be complicated, this would favor surveys 
like HETDEX and BOSS versus WiggleZ. 

Finally, the choice of a different split in k space mildly 
affects the resulting Figure of Merit and the correlation 
between e and v, as long as the achieved precisions in 
the two bins are not very different (in which case the 
uncertainties become larger). One could choose k mc d in 
different ways, for example choosing the one that mini- 
mizes the FoM or the correlation between the two param- 
eters. However, this optimization works for one survey 
at a time. The prescription used here to determine fc me d 
allows for a simple and fair comparison of multiple sur- 
veys. 

For the same reason, we chose to focus in this paper on 
the minimal number of bins in wavenumber and in red- 
shift that allow both parameters to be measured. The 
combined power of surveys at different redshifts to dis- 
criminate against GR+ACDM can be easily obtained by 
adding their Fisher Matrices. There are obvious advan- 
tages in having a broad lever arm in both the radial (z) 
and transverse (fc) direction; if the concordance model 
were falsified, breaking down the information might re- 
veal deviations in only one bin, and cast light on the 
nature of such deviation. However, looking at probes 
at multiple redshifts and/or several bins in k is survey- 
specific and is beyond the scope of this paper. 

V. CONCLUSIONS 

We have shown how redshift galaxy surveys can 
be used to test the consistency of the concordance 
GR+ACDM model with the observations (US]. 

This approach is motivated by the idea that the sim- 
plest model should be abandoned only if evidence of its 
failure is found. This framework allows one to search 
for deviations from the concordance model behavior in 
a completely model-independent fashion. The two fea- 
tures of the concordance model we explored are the fixed 
relationship between the expansion history and the log- 
arithmic derivative of the growth, and the fact that the 
latter is independent of scale. Each of these relation- 



ships can be rewritten in terms of a parameter whose 
value is known in the concordance model. We have de- 
fined these parameters as e and v, and used a Fisher 
Matrix formalism to derive a Figure of Merit for three 
galaxy surveys, WiggleZ, BOSS and HETDEX, in the 
(e, v) plane. For this particular choice of observables and 
techniques {e.g., the use the bispectrum to determine the 
bias), and assuming that any scale dependence of bias can 
be removed, high- redshift, low-bias galaxies are found to 
be more powerful in detecting the presence of deviations 
from the concordance model. However, if such deviations 
were found, characterizing them would require a compre- 
hensive mapping of the expansion history and growth of 
structure on multiple redshifts and scales. 

The framework we developed for testing the 
GR+ACDM model can be easily applied to other mea- 
surements of e and v; many of current constraints on 
the growth of structure |34l50| could be reinterpreted 
using this formalism. We believe that this method is 
useful for two reasons. On the one hand, while it is not 
conceptually different from fitting the parameters of the 
GR+ACDM model, it enables straightforward, separate 
tests of two of its properties. On the other, the Figure of 
Merit we have introduced allows one to compare differ- 
ent surveys on an equal footing, and to easily predict the 
combined power of different surveys in detecting devia- 
tions from the concordance model, which is a new feature 
of this approach. 

The method can also be extended to different ob- 
servables. For example, in non-standard cosmological 
models, the scale dependence of the growth factor D(k, a) 
and its logarithmic derivative / can be completely dif- 
ferent. As a result, a simultaneous measurement of these 
two quantities would be another model-independent test 
of the concordance model. Galaxy redshift surveys can 
in principle measure both, since the redshift distortion is 
sensitive to /, but the power spectrum itself is sensitive 
to the growth factor D; the latter can also be measured 
via weak gravitational lensing [5lH54j |. and galaxy 
clusters counts [55j . We defer this extension to another 
paper [56j . 
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Hill, Karl Gebhardt, Eiichiro Komatsu, and David 
Spergel for useful comments, and Donghui Jeong and 
Masatoshi Shoji for enlightening discussions. We also 
thank the referee, Eric Linder, for many suggestions 
which helped us improve the manuscript. 
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